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We calculate the energy of a Yang-Mills vortex as function of its magnetic flux or, else, 
\ of the Wilson loop surrounding the vortex center. The calculation is performed in the 

1-loop approximation. A parallel with a potential as function of the Polyakov line at 
nonzero temperatures is drawn. We find that quantized Z(2) vortices are dynamically 
Q_i' preferred though vortices with arbitrary fluxes cannot be ruled out. 



The hypothesis that Z(N C ) vortices are responsible for confinement has several attractive 
features. First, one immediately gets the area law for Wilson loops in the fundamental 
representation of the SU(N C ) gauge group, even under the simplest assumption that vortices 
are non-interacting and therefore the number of vortices piercing a given loop is Poisson- 
distributed Jl|. Furthermore, Wilson loops in any representation transforming under the 
group center Z(N C ) will have the same string tension as in the fundamental representation, 
while the 'center-blind' representations will have, asymptotically, zero string tension. This 
is what is expected from screening by gluons. The temporal approximate Casimir scaling of 
the string tensions may be probably explained by the finite sizes of the vortex cores 0. 

Second, heating the d — 3 + 1 Yang-Mills system towards the high-temperature d = 3 
case and then heating d = 2 + 1 system eventually to the d = 2 case, one does not observe 
abrupt changes in the spatial string tensions. This continuity as one goes from four to two 
dimensions favours d = 2 objects as being basic for confinement in all the above cases []. In 
three dimensions vortices form closed lines; in four dimensions they form closed surfaces. 



In d — 2 pure Yang-Mills theory the confinement is trivial in axial gauges and may seem to have no 
relation to vortices. However, it is known that at least on compact d = 2 manifolds the theory is exactly 
equivalent to a sum over vortices, see a recent paper B and references therein. 



For vortices to be physical objects and not artifacts of a regularization, their cores should 
be finite in physical units, i.e. to be of the order of 1/A ~ M _1 exp(3 • 87r 2 /HN c g 2 ) in d = 4 
and of the order of l/g 2 . in d = 3. Correspondingly, the energy of a vortex per unit length 
should be of the order of g\ in d = 3, and the energy per unit surface should be of the order 
of A 2 in d = 4. First indications that in might indeed be the case came recently from lattice 
studies using the smoothing procedure ||. 

In order to reveal 'thick' vortices theoretically, one has first of all to integrate out the 
high-momentum components of the Yang-Mills field. The resulting effective action may then 
have a stable saddle point of a vortex type, with the core size remaining finite as one takes 
the UV cutoff to infinity. When and if it happens, one can speak of vortices as physical 
entities, and further on to try to build a theory of those extended objects in the vacuum. 

The main dynamical variable of a vortex is the azimuthal component of the Yang-Mills 
field A^(p) = e a/ 3n a Ap, where n a is a unit vector in the plane transverse to the vortex. 
One can always choose a gauge where is independent of the azimuth angle 0. Generally 
speaking, it implies that the radial component A p (p) ^ 0, however, we shall neglect this 
component as it can be always reconstructed from gauge invariance by replacing d p — ► 
d p 5 ab + f acb A c p . A circle Wilson loop lying in the transverse plane and surrounding the 
vortex center is then 

Wi(p) = ^TrPexpifpA^ = cos[7r/i(p)], p{p) = pj A%{p)A%{p), (1) 

taking for simplicity the fundamental representation of the SU(2) gauge group. For an 
arbitrary representation of SU(2), labelled by spin J, one has 

2 J + 1 sm[7ip{p)\ 

If p(p) — > integer at large distances p from the vortex core, the Wilson loop Wj(p) — > (— 1) 2J . 
This is the definition of the Z(2) vortex. We shall see that integer values of p(oo) are 
dynamically preferred. 

Neglecting all components of the Yang-Mills field except the essential azimuthal one, the 
classical action of the vortex becomes 



d d X^f- = ld d ~ 2 X, [d 2 X L [ -^f- = {d d ~ 2 X\\ 27T / 
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dp p 



(3) 



To get the effective action for p(p) we integrate over quantum fluctuations about a given 
background field A^(p)p considered to be a slowly varying field with momenta up to certain 
k min . Accordingly, quantum fluctuations have momenta from k min up to the UV cutoff k max . 



Writing A, 



Af, + a M where A p 



Sf^A^p) we expand F 2 U (A + a) in the fluctuation field 



a M up to the second order appropriate for the 1-loop calculation. The term linear in 
vanishes due to the orthogonality of high and low momenta. The quadratic form for a M 
the standard 



is 



W«l = -[D 2 (A)] a \ u - 2f acb F c (A) 



(4) 
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if one imposes the background Lorentz gauge condition, 

Df{A)a b ^ = 0, Df{A) = d,6 ab + f^A^. (5) 

The effective action is 

S eff [A] = ilndet(^) -lndet(-TJj). (6) 

In the presence of dynamical fermions in the fundamental representation the Dirac de- 
terminant should be added: 

- In det( V m7 „) = -\ In det (vj - \ > V m = «m ~ ^m^- ( 7 ) 

The effective action may be expanded in powers of the (covariant) derivatives of the 
background field A^. We are now interested in the first nontrivial term of this expansion, 
namely in the zero-derivative term, that is in the effective potential as function of the flux 
pA§. The next term with two derivatives will renormalize eq. (|3|), and there will be, generally 
speaking, further terms. In the zero-derivative order F tlu (A) = 0, therefore det(W Atl/ ) = 
d-det(— D 2 ), where d is the full space dimension. With A^ being the only nonzero component 
of the background field, it is natural to write down D 2 in the cylindric coordinates: 

The effective action ([]) for slowly varying pA^ is 



SefAM = - \± - lj Sp exp(sD 2 ) = Jd d - 2 x {l Jd 2 x ± V^W, /i = P^A^, 

(9) 

where we have introduced the 'vortex potential energy' in d dimensions V^(p). Here Sp 
denotes the functional and the colour traces, and integration over the proper time s has 
been introduced. The limits s min ~ l/k"^^ and s max ~ l/k^ n/in are gauge- invariant UV and 
IR cutoffs, respectively. In saturating the functional trace one can use any complete set of 
functions. Naturally, one uses the plane wave basis for the longitudinal components x^.-.Xd, 
and the polar basis exp(zm0) fk,m(p) f° r the transverse ones, where the functions fk,m(p) 
must form a complete set, see below. In what follows we shall consider only the 577(2) 
case. The transverse part of the covariant Laplacian (§), after acting on the polar harmonics 
exp(im0), m = integer, has three eigenvalues as a 3 x 3 colour matrix: 

1 d d if m2 ' 

-frpf> — H ( m +/ i ) 2 ' ( 10 ) 

pdp dp p 2 | (m _^ )2 _ 

The first one cancels out when one subtracts the free determinant (without the external 
field); the last two are actually coinciding, given that m ranges from minus to plus infinity. 
The differential operator (eq. fllOD , last line) has eigenfunctions J±( m -n)(kp) with eigenvalues 
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—k 2 . The index of Bessel functions must be non- negative to ensure regularity at the origin. 
We thus choose the functions -F m fc (0, p) = exp(im0) J\m-^\{kp) as a complete functional basis 
in the transverse plane: 



1 

2^ 



E 



m=— oo 



oo I 

dk k F m ,*(0, P )F* >fe (0', p') = -8{p - p') 5(<p - 0') l mod 2 , 



These functions satisfy also the ortho-normalization condition: 

2vr Jo ^ Jo dp 9 Fm ' fc ^' P) = k 6< " k ~ k "> 5mm '- 



(11) 



(12) 



Eqs.(|TT], 12) can be checked using the formula 6.541.1 from Gradshteyn and Rhyzhik ||. 
Eq. (0) can be thus rewritten as 



V {d \p) = -(d-2) 



dp 3 ...dp d 



%ax ds f 00 



S JO 



(27t) d - 2 

•gZ E [ J \m-^\( k P) - J \m\( k P) 



k dk exp[— s(k 2 + pi + ... + p^)] 



(13) 



We see that the dependence on the flux p enters only through the indices of the Bessel 
functions and that the potential is explicitly periodic in p with period 1. Integration over 
momenta can be easily performed using eq. 6.633.2 from |J, yielding 



Smax dS f 1 



V^(p) = -(d-2) ™ 
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exp 



m=— oo 



2s 



I\m\ I — 



2s 



(14) 



We next sum over m using the integral representation for modified Bessel functions, eq. 
8.431.5 of ||. To write explicit formulae we imply that p e (0, 1). If p is outside this 
interval V i - d \p) should be continued by periodicity. We get 



V (d \p) 



d-2 sin(vr/i) f°° dx cosh (§ ~~ p) x 



(27T P 



2^12 



7T 



cosh | 



"dt^ _1 exp [-t(coshx + 1)] , (15) 



where we have introduced a new variable t = p 2 /(2s) instead of s. Correspondingly, the 
integration limits become t min = p 2 /2s max and t max = p 2 /2s min . It should be stressed that 
the above expression is finite in the limit when one removes both the UV cutoff (t max — > oo) 
and the IR cutoff (t m ir, 
we finally get: 



V {d \p) 



0). In this limit both integrations in (15) become elementary, and 

d-2 sinM r (f ) r (f + p) r (f + 1 - p) 

(np 



d 

2^ 2 



p 4 12 
1 l tan(7Tju) 



p 6 96 



\s p{l - p 2 ){2 - p) 
l-4p 2 )(3-2/i) 



mod 1 



r(d + i) 

for d — 4, 



(16) 



mod 1 



for d = 3. 
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Energy as function of vortex flux 
in three and four dimensions 
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Figure 1: The vortex potential, eq. (|T6|). 

At d = 2 the vortex potential is identically zero as there are no transverse gluons in two 
dimensions. Though the potential is given just by a polynomial (in case d = 4), it is actually 
a periodic function of \i with a unit period, as seen explicitly from eq. (|T3|). At integer \i the 
potential is zero but has a jump in the first derivative; it is depicted in Fig.l. 

The problem we have solved has certain similarity with that of the potential energy as 
function of the d-th component of the Yang-Mills field A° =(J , in case of nonzero temperatures. 
In that problem one integrates over all Matsubara frequences in the background of a static 
A a dl where d = 4 for 3 + 1 dimensions, and d = 3 for 2 + 1 dimensions. The two problems 
are ideologically similar, only the topology is different: in the nonzero-temperature case the 
topology is that of a cylinder x with the compact dimension in the 'temperature' 

direction, while in the problem considered here the topology is that of a plane with a deleted 
point at the origin, \ {0} x R( d ~ 2 \ The role of Matsubara frequencies is played by 
the polar harmonics exp(im0). In the first case one finds the potential energy as function 
of the Polyakov line winding in the compact dimension; in the second case one finds the 
potential energy as function of the Wilson loop winding around the point at the origin. In 
both cases the potential energy is a periodic function, however the calculations are easier in 
the nonzero-temperature case as one can evaluate the functional trace (|9]) in the plane-wave 
basis. A simple calculation of eq. @ gives the following form of the nonzero-temperature 
potential energy as function of v = J A d A d /(2%T), valid for any number of dimensions: 
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(17, 

It is explicitly periodic in z/. At d = 4 the integration can be performed with the help of eq. 
3.533.3 [H, and we get the well-known result ||: 



\2 



12tt 2 



= ^K) 2 + 0(A 3 ). (18) 

mod 1 



At d = 3 the integral ( |TT| ) cannot be expressed through elementary functions; however it can 
be compactly written as 

2T 3 . o r°° , x 2 cosh x 



21/ Z 100 
\z( 3 )(zy) = sin 2 7T^ / dx 

IT JO 



sinhx(sinh x + sin ~kv) 



= In 1 + 0(z/) = £(^) 2 In ^ + 0(A 4 ). (19) 

To our best knowledge, this is a new result. 

It is worth mentioning that for d = 3 the Debye mass is infrared-divergent, hence the 
non-analiticity in the v 2 term. In the case of the vortex potential ([if]) the infrared divergency 
is even stronger: as a result the expansion of the potential starts from the non-analytic 
term both for d = 3 and 4. It is interesting that in all cases above the non-analytic terms are 
due to the contribution of zero Matsubara frequencies (in case of the vortex "zero frequency" 
means the azimuthal harmonic with m — 0). Indeed, the contributions of zero frequencies 
to the 'temperature' potentials ( |T8"D and flT9| ) are 

V^%{v) = - ^J] 2u 3 + UV divergence, V^%(v) = ttTV In \ + UV divergence. 

(20) 

The m = contribution to the vortex potential ( fL6|) is 

Vl 4 2 (^) = \zr^{/i - /i 3 ) + UV divergence. (21) 

These are exactly the non-analytic terms in all the cases. The UV divergences are cancelled 
by contributions from nonzero frequencies. 

If there are dynamical fermions in the problem, they should be treated differently in the 
nonzero-temperature case and in the vortex one. In the first case fermions are antiperiodic in 
the 'temperature' direction, hence one has to sum over half-integer Matsubara frequencies. 
In the latter case fermion wave functions are periodic functions of the azimuthal angle 0. The 
resulting fermion contribution to the vortex potential (in the fundamental representation) 
is obtained from eq. fll6D by substituting \i — > fi/2 and multiplying the result by —2. The 
fermion contribution is periodic in // with period 2, and not 1 as in the boson case. It can 
be easily checked that in the supersymmetric case, with Majorana gluinos belonging to the 
adjoint representation, the fermion contribution cancels exactly with the boson one, so that 
the vortex potential is zero. 
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Returning to the bosonic contribution to the effective vortex potential ( |I6|) in the pure 
glue case, one concludes that integer values of /i(p) are a must at p — > 0, otherwise the 
integral over p, 



diverges. Integer values of p(p) at p — ► oo are, clearly, energetically favourable though the 
energy fl22|) remains finite at noninteger values of p(oo) as well because of the convergence of 
the integration over p at large p. It means that the quantized Z(2) vortices are dynamically 
preferred but noninteger fluxes are not altogether ruled out by the energetics. This should 
be contrasted to the case of nonzero temperatures where quantized values of Ad at spatial 
infinity are necessary to make the energy finite. 

I am grateful to Victor Petrov, Maxim Chernodub and Konstantin Zarembo for useful 
discussions. 
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